On the study of jamming percolation 
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We investigate kinetically constrained models of glassy transitions, and determine which model 
characteristics are crucial in allowing a rigorous proof that such models have discontinuous transi- 
tions with faster than power law diverging length and time scales. The models we investigate have 
constraints similar to that of the knights model, introduced by Toninelli, Biroli, and Fisher (TBF) , 
but differing neighbor relations. We find that such knights-like models, otherwise known as models 
of jamming percolation, need a "No Parallel Crossing" rule for the TBF proof of a glassy transition 
to be valid. Furthermore, most knight-like models fail a "No Perpendicular Crossing" requirement, 
and thus need modification to be made rigorous. We also show how the "No Parallel Crossing" 
requirement can be used to evaluate the provable glassiness of other correlated percolation models, 
by looking at models with more stable directions than the knights model. Finally, we show that the 
TBF proof does not generalize in any straightforward fashion for three-dimensional versions of the 
knights-like models. 



I. INTRODUCTION 

The puzzle of glass-forming systems has remained suf- 
ficiently elusive over the years such that even the puz- 
zle pieces themselves have changed shape. For exam- 
ple, the puzzle piece of the lack of a growing lengthscale 
may now have to be modified since a growing lengthscale 
can perhaps be extracted from a higher-order correlation 
function A less recent change in puzzle pieces is the 
distinction by Angell between fragile and strong glasses, 
where certain effects are more dramatic in fragile glasses 
than in strong ones 2]. One piece of the puzzle that 
has remained constant over the years, however, is the 
dramatic slowing down of the dynamics of the particles 
near the glass transition. More precisely, a supercooled 
liquid's viscosity can increase by fourteen orders of mag- 
nitude as the temperature is decreased near a "working" 
definition of the glass temperature 0]. 

Two main phenomenological models for this dynam- 
ical slowing down have emerged over the years: mode- 
coupling theory and kinetically constrained models. We 
will not focus on mode-coupling theory here and sim- 
ply refer the reader to several excellent reviews [!, 0, H[ • 
As for the second phenomenological approach, kinetically 
constrained models, the goal is to understand whether 
glassy dynamics can be understood as arising from steric 
constraints on the particles alone Q . One of the simplest 
such examples is the Kob- Andersen model [7] . It is mo- 
tivated by the caging of particles, ultimately observed in 
larger scale systems such as colloidal glasses 18, 9]. The 
Kob-Andersen model is a hard-core lattice gas model, 
but with the constraint that a particle can only hop to 
an adjacent site if and only if it has less than a certain 
number of neighbors, m, both before and after the move. 
Early simulations of the Kob-Andersen model on the hy- 
percubic lattice for relevant values of m appeared to find 



* mjeng@physics.syr.edu 
tjschwarz@physics.syr.edu 



a dynamical phase transition at a nontrivial critical den- 
sity 0]. However, subsequent mathematically rigorous 
results found that the phase transition does not occur 
until the fully packed state [l(J, EH ■ This corresponds to 
a zero-temperature glass transition. TBF proved this by 
showing that at any monomer density, there were mobile 
cores that could diffuse at sufficiently long time scales. 

While the hypercubic version of the Kob-Andersen 
model does not exhibit a finite-temperature glass transi- 
tion, the mean field version does [l(J, [TTj] . Since it is still 
up for debate whether or not mean field is relevant for 
physical systems, one can ask whether or not there exists 
another finite-dimensional, kinetically constrained model 
that exhibits a finite temperature glass transition. While 
the Kob-Andersen (and the Fredrickson- Andersen (I2I ]) 
models, are elegant in their simplicity, there are in- 
deed two more involved kinetically constrained models 
in two dimensions that can be proven to exhibit a finite- 
temperature glass transition [3, EE EE]- These two 
models have been dubbed the spiral model [I?) and the 
sandwich model [l8|. Both models exhibit an unusual 
phase transition in that the fraction of frozen particles 
jumps discontinuously at the transition, typical of a first- 
order phase transition. However, as the transition is ap- 
proached from below, there exists a crossover lengthscale 
that diverges faster than a power law in T — T g . The 
crossover length, T, distinguishes between squares of size 
L << r, which are likely to contain a frozen cluster, 
and squares of size L » T, for which the probability 
of containing a frozen cluster is exponentially unlikely. 
Given this combination of a discontinuity in the fraction 
of frozen particles, and a faster than power law diverg- 
ing length scale, the transition has unique characteris- 
tics pjf. 

Models such as the spiral and sandwich models are 
proof in principle that further exploration of kinetically 
constrained models in finite-dimensions may be fruitful, 
in particular, because should an ideal glass transition ex- 
ist, it may indeed be of unusual character. More specif- 
ically, the Edwards-Anderson order parameter should 
be discontinuous at the transition, yet accompanied by 



2 



rapidly diverging time scales [H, H(| ■ 

Kinetically constrained models can be related to mod- 
els of correlated percolation, which are models in which 
each site is initially occupied with an independent prob- 
ability p, as in normal percolation, but then correlations 
are induced through some culling rule — that is, sites can 
only be occupied if certain conditions on their neighbor- 
ing sites are met. The two types of models are related, 
with the immobile particles of the kinetically constrained 
models corresponding to the stable occupied particles of 
the correlated percolation models. The simplest model 
of correlated percolation is fc-core percolation, in which 
an occupied site must have at least k occupied neigh- 
bors [lo, [Sj]. Occupied sites that do not satisfy this 
stability requirement are removed, and this condition is 
applied repeatedly, until all remaining sites are stable. 
This model can be mapped to the Fredrickson- Andersen 
model [EUl. 

One might think that kinetically constrained mod- 
els and correlated percolation systems would be easy to 
numerically simulate, and that construction of rigorous 
proofs would be nothing more than an interesting prob- 
lem for mathematicians. However, while simulating these 
systems is easy, extracting their properties in the infinite 
system size limit is much harder. Just as with the Kob- 
Andersen model, initial numerical simulations of fc-core 
percolation for certain values of k found evidence of first- 
order phase transitions at nontrivial critical densities, 
and of second-order phase transitions in a different uni- 
versality class than normal percolation [23|, [24|, [2|| . How- 
ever, subsequent mathematically rigorous analyses found 
that the critical point is p c — 1 for those k [2y, [U [H[ • 
Because the critical point only approaches unity very 
slowly in the limit of infinite system size, the simula- 
tions on finite-size systems were misled as to the location 
of the true critical point, thereby highlighting the impor- 
tance of rigorous results for these models as well. For a 
review of fc-core percolation, see Ref. [29[ . 

There is another finite-dimensional system exhibiting 
an unusual transition, at least numerically. It is the 
jamming transition in repulsive soft spheres. Numeri- 
cal simulations of soft spheres show a critical point at 
which the average coordination number jumps discon- 
tinuously to a universal, isostatic value. But quantities 
such as the shear modulus and the deviation of the av- 
erage coordination number from its isostatic value show 
a nontrivial power law behavior in the vicinity of the 
critical point [3ll . I32I [33|. Recent experiments on two- 
dimensional photoclastic beads support this notion of a 
mixed transition [34J. 

Interestingly, it has been conjectured that the physics 
of granular systems, colloidal systems, and glassy systems 
are of a similar character [35[ . The mean field results of 
not one but several correlation percolation models corre- 
sponding to glassy physics support this notion quantita- 
tively. Furthermore, experimental evidence of caging in 
another two-dimensional granular system also supports 
this notion [36] . The question of finite-dimensional glassy 



models being quantitatively similar to the repulsive soft 
sphere system is still being investigated. Certainly the 
spiral and sandwich models show that, qualitatively, one 
can have a glassy system exhibiting an unusual finite- 
dimensional transition. However, they do not appear to 
be in the same class as the jamming system, since the 
order parameter exponent is unity just above the transi- 
tion in the jamming percolation models, but is one- half 
in the jamming system. 

To explore the possible link between jamming and 
glassy systems in terms of a finite-temperature glass 
transition, TBF initially introduced the knights model, 
a model of correlated percolation similar to the spiral 
model [3, and called it a model of jamming percola- 
tion. In fact, the spiral, knights, and sandwich models 
are all models of jamming percolation. In this paper, we 
expand on our earlier work (Ref. [Hj]), in which we in- 
troduced the sandwich model, by presenting the details 
of the proof of an unusual transition in this model. This 
proof is based on modifying the proof developed by TBF 
in Refs. [3, HE EE] for the spiral model (although orig- 
inally misapplied to the knights model). We also intro- 
duce further generalizations of models of jamming perco- 
lation to demonstrate that the phenomenon of a finite- 
dimensional transition is indeed somewhat generic. In 
doing so, we show that the TBF proof only gives a rigor- 
ous derivation of these novel properties if a "No parallel 
crossing" rule holds. This rule says that two similarly- 
oriented directed percolation chains cannot cross without 
having sites in common. The effect of this rule is that 
one directed percolation-like process cannot be used to lo- 
cally stabilize the other. For models such as the sandwich 
model, which satisfy this "No parallel crossing" rule, but 
fail a "No perpendicular crossing" rule, the TBF proof 
works only with some modification. 

The methods described in this paper can be used to un- 
derstand for which correlated percolation models a proof 
along the lines of the TBF proof can be used to show 
a glassy transition, and for which models they cannot. 
Given that there are two very detailed papers on the 
TBF proof [TBI, H6| . we will refer to them quite often, 
as opposed to making this paper self-contained. Finally, 
we will discuss connections between jamming percolation 
and force-balance percolation, another correlated perco- 
lation model inspired by granular systems, where numer- 
ical evidence points toward an unusual transition in finite 
dimensions [37l. [38j. 



II. THE CLASS OF MODELS 

We consider a class of models that generalizes the 
knights model, the earliest of the jamming percolation 
models. The class of models is defined on the two- 
dimensional square lattice. Initially, each site is occupied 
with an independent probability p. Each site has four 
neighboring sets, and each set contains two sites. The 
four sets are labelled as the northeast, northwest, south- 
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FIG. 1: The knights model neighbors. 
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FIG. 2: The sandwich model neighbors 



east, and southwest neighboring sets. (The sites in those 
sets only lie in precisely those compass directions for the 
knights model, but we continue to label the four sets in 
this manner for all our correlated percolation models un- 
til section[lV]) To be stable, an occupied site must either 
have (1) at least one northeast neighbor and at least one 
southwest neighbor, or (2) at least one northwest neigh- 
bor and at least one southeast neighbor. All other oc- 
cupied sites are unstable, and are vacated. This culling 
process is then repeatedly applied — sites that were pre- 
viously stable may become unstable by earlier cullings — 
until all remaining sites are stable. The neighboring sets 
of the original knights model introduced by TBF [lj] are 
shown in Fig. [T] Fig. [3 shows the neighboring sets in the 
sandwich model, which we introduced in [l8[ , and Fig. [3] 
shows the neighboring sets in the spiral model, which 
TBF introduced in 

In these models, stable sites must either be part of a 
chain running from the northeast to the southwest, or 
part of a chain running from the northwest to the south- 
east. In the final configuration, all sites must be stable, so 
any chain must either continue forever (to the boundary) , 
or terminate in a chain of the other type (this latter case 
is called a "T-junction"). Any sites left after the culling 
procedure must thus be connected by an infinite series 
of chains (or, in a finite system, connected by a series of 
chains to the distant boundary), so asking for the criti- 
cal probability at which an infinite cluster first appears is 
the same as asking for the minimum probability at which 
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FIG. 3: The spiral model neighbors 
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FIG. 4: A sublattice in the sandwich model. 



some sites remain unculled in the infinite size limit. 

If sites could only be stable by having northeast 
and southwest neighbors, then every stable site would 
be part of a northeast-southwest chain on a particular 
sublattice — the sublattice extending to the northeast of 
the site for the sandwich model is shown in figure 2] 
Chains on this sublattice are isomorphic to infinite chains 
in directed percolation, which has a well-studied second- 
order phase transition 39]. However, in these models 
there is an additional mechanism by which sites can be 
made stable — that is, by having northwest and southeast 
neighbors. Adding an extra way for a sites to be stable 
can only possibly depress the critical probability, so we 
immediately see that for these models p c < p^ p . 

It will turn out to be useful to divide these models into 
classes based on two properties. We define a model as 
having "No Parallel Crossing" property if whenever two 
northeast-southwest chains (or two northwest-southeast 
chains) intersect, they must have sites in common — we 
abbreviate this a "Property A." And we define a model as 
having a "No Perpendicular Crossing" property, or Prop- 
erty B, if whenever a northeast-southwest and northwest- 
southeast chain intersect, they must have sites in com- 
mon. Table [T] shows which properties each of the three 
models possesses: 



Model 


Property A 


Property B 


knights 


No 


No 


sandwich 


Yes 


No 


spiral 


Yes 


Yes 
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(a) 



(b) 



(c) 



FIG. 5: (a) Failure of Property A ("No Parallel Crossings") 
in the knights model, (b) Failure of Property B ("No Per- 
pendicular Crossings") in the knights model, (c) Failure of 
Property B in the sandwich model. 



Examples of where the properties fail for each model are 
shown in figure [5] 

Our analysis here is based to a large extent on the 
claimed TBF proof of a glassy transition for the knights 
model. Their proof consisted of two parts. First, they 
claimed to show that the critical point of the knights 
model is exactly the same as that for directed percola- 
tion. Second, once this was done, they were able to use 
well-known results on directed percolation (assuming a 
well-tested conjecture about anisotropic rescaling in di- 
rected percolation) to show that this model has a glassy 
transition — specifically, they were able to find structures 
with a finite density at the critical point of directed per- 
colation, and to show that just below this critical point, 
the crossover length and culling times di verg ed. A short 
version of their proof appeared in Re fs. Jl4l| . with more 
detailed explanations in Ref. [l5| and (lq |. 

Our analysis of this more general class of models shows 
that the TBF proof that p c — p^ F is only valid for mod- 
els satisfying property A; so it works for the sandwich 
and spiral models, but fails for the knights model. The 
second part of their proof, showing a glassy transition, 
implicitly assumes property B. The spiral model exhibits 
property B and hence the TBF method of proof carries 
through [40(. However, we show that the proof can be 
modified to work for models that fail to have property B. 
The spiral and the sandwich models thus have provably 
glassy transitions, while the knights model does not. 



III. IDENTIFYING THE CRITICAL POINT 

We sketch the TBF proof that plights = p DP ^ which 
will let us understand why property A is sufficient, and 
most likely necessary to the result. The key to the TBF 
proof is to show that voids (clusters of empty sites) of 
particular shapes have a finite probability of growing for- 
ever. For example, for the diamond-shaped void in the 
sandwich model, shown in figure [BJ if the key site la- 
belled y is vacant, it will trigger the removal of all the 
sites marked with stars, increasing the void size by one. 
(The corresponding void for the knights model appears 
as figure lc of [14].) If this process repeats forever, with 
such key sites repeatedly removed, this void will grow to 
infinity. 
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FIG. 6: A void in the sandwich model. If the site y is unsta- 
ble, its removal will cause the culling of all the sites marked 
with stars. The corresponding void for the spiral model is a 
diamond, regular in shape. 



The TBF proof for the knights model is based on the 
claim that the key sites of the knights model, located at 
the corners of octagonal voids of size L, can only be sta- 
ble if part of a directed percolation chain of O(L). Once 
this claim is granted, the rest of the proof is straightfor- 
ward. For p < p^ p , such long chains are exponentially 
suppressed, and thus for a large void, the vertices at the 
corners of the void are exponentially likely to be culled. 
Summing up the relevant probabilities, this results in a 
finite probability that the void will grow to infinity. For 
an infinite lattice, it is thus certain that there will be 
at least one void that grows to infinity, showing that all 
sites are culled below pj? p . Since plights < p® p i this is 
supposed to show that p^ m s hts _ ^dp 

The claim that sites at the corners of voids can only 
be stable if part of a long chain of O(L) is only valid 
for models with property A. A counterexample to this 
claim for the knights model can be seen in Fig. [7J This 
counterexample also shows why property A is necessary 
and sufficient for this claim to be valid. To be stable, the 
site must be part of a northeast-directed chain; pick the 
lowest northeast-directed chain coming out of the corner 
site y. If that chain stops before reaching length O(L), 
it must terminate in a northwest-southeast chain. Since 
the northeast chain is not as long as a wall of the void, 
the new southeast-directed chain will eventually hit the 
void (thus resulting in the culling of all chains, and the 
corner site y), unless it hits a T-junction. That new 
T-junction results a second northeast-southwest chain, 
which will eventually reach the first northeast chain, as 
in figure [71 For models with property A, the two chains 
will intersect, contradicting the original assumption that 
we chose the lowest northeast-directed chain out of y. 
Thus, by contradiction, for models with property A, the 
first northeast chain must be O(L) for y to be stable, and 
it indeed follows that p c = pj? p . 

What about for models such as the knights model, 
that lack property A? Is it possible that despite this 
counterexample to the claim, pkmg hts j g ac tually equal 
to pj? p , for some other reason? While we do not have 
a mathematically rigorous proof that pkmg hts _^ ^dp ^ we 
present an argument that the two are almost certainly 
unequal. We present our arguments in the context of the 
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FIG. 7: A counterexample to the claim in Refs. [bH . HH Il6| ] 
that for y to be stable it must be part of a long, uninterrupted 
chain to the northeast. Large solid circles are occupied sites, 
and large empty circles are vacant sites. All sites in the void 
are vacant. All other sites can be either occupied or unoccu- 
pied. Arrows are drawn from each occupied site to indicate 
the neighboring sites that make it stable. 
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FIG. 8: A substructure that depresses the critical point of the 
knights model. Numbers indicate sites in the same sublattice. 



knights model, but they generalize to other models that 
lack property A. 

Consider the substructure in Fig. [5] All sites in it are 
stable under the knights model culling rules, except for 
the two sites at its ends, and those sites will become sta- 
ble if the substructure is attached between two northeast- 
southwest chains. Furthermore, there is no northeast- 
southwest chain internal to the substucture connecting 
the two ends. The substructure is internally stabilized 
by northwest-southeast links. 

This means that the substructure can act as a "rest 
stop." Northeast-southwest chains can have breaks in 
their paths, connected by this substructure. Just below 
p^ F , long northeast-southwest chains almost form an in- 
finite structure. They are almost linked, so a few extra 
connections, through these substructures, should create 
an infinite cluster even below p® p . So we expect that 
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FIG. 9: A configuration in "jumping directed percolatiom 
that receives an additional connection with probability s. 



knights DP 
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We can make the argument more formal by consider- 
ing the following modification of the directed percolation 
problem, which we call jumping directed percolation. As 
with normal directed percolation, we occupy sites on the 
square lattice with probability p, and connect each site 
with directed bonds to its neighbors to the north and 
east. However, now we define an additional way for sites 
to be connected. We divide the lattice into blocks of 
size 9x9, and for each block, if the two hollow squares of 
sites shown in Fig. [5] are occupied, we with probability 
s connect the two hollow squares with a directed bond 
from the southwest square to the northeast square. The 
critical point of this model is a function of s: p^ ump (s). 

By repeating Fig. [8] three times, to create diamonds 
connecting sublattice #1 to #2 to #3 and then back 
to #1, we obtain a structure that links two separated 
diamonds on sublattice #1, through sites in the other two 
sublattices. So if we restrict ourselves to looking at sites 
in sublattice #1, and the directed percolation structures 
on that sublattice, sites that appear disconnected may be 
connected by these sites in sublattices #2 and #3. The 
structure obtained by repeating Fig.[8]three times has 24 
sites in sublattices #2 and #3, each of which is occupied 
with probability p, and has 24 sites in sublattice #1. The 
sites in sublattice #1 in this repeated structure map onto 
the occupied sites in Fig. [9] So with s = p 24 , s gives the 
probability of having appropriate "hidden" occupied sites 
in sublattices #2 and #3 that connect and make stable 
the two hollow squares. Infinite chains in the jumping 
directed percolation model are infinite stable clusters in 
the knights model, and thus ^nights < j0 Jum P ^24) < p DP 

However, since jumping directed percolation is just di- 
rected percolation with an extra connection process, it is 
reasonable to expect that p^ ump (s) < p^ F for all s > 0, 
implying pj™ights < p® F . While this argument is not 
mathematically rigorous, it is strongly suggestive, partic- 
ularly when we recall previous results on enhancements in 
percolation by Aizenman and Grimmett (4lj . Their work 
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showed that if percolation on a lattice was "enhanced" by 
adding, for specified subconfigurations, extra connections 
or occupations with probability s, this would strictly de- 
crease the critical probability, for any s > 0, so long as the 
enhancement was essential. Essential enhancements were 
defined as those such that a single enhancement could 
create a doubly-infinite path where none existed before. 
See Ref. [4l| for a more rigorous and precise statement of 
the results on enhancements, and Ref. [42j for a general 
discussion of enhancements. The results of Ref. [4jJ were 
obtained for undirected percolation, so are not directly 
relevant for the jumping directed percolation model con- 
sidered here, but they are analogous enough to strongly 
suggest that p J c ump (s) < p® p for all s > 0. It is diffi- 
cult to see how such adding such a new route for paths 
to infinity could leave the critical probability completely 
unchanged. 



IV. PROPERTY B 

For models satisfying property A, we have p c = p^ p , 
but we still need to check that the TBF proof that 
the transition is glassy (discontinuous with a diverging 
crossover length) is valid. 

The TBF proof of a glassy transition implicitly as- 
sumes that the knights model has property B. For ex- 
ample, to show discontinuity, they construct a configura- 
tion that has a finite density at pj? p — see Fig. 2b of [13] • 
This figure, and others like it, implicitly assume prop- 
erty B, because they are based on drawing overlapping 
rectangles in independent directions, and assuming that 
if paths in these intersecting rectangles cross, they must 
stabilize each other (form a T-junction). The resulting 
frozen structure is shown on the left side of figure 1101 
However, if a model lacks property B, the paths can cross 
without stabilizing each other. 

The knights model does not satisfy Property A, so 
whether or not it satisfies Property B is a moot point. 
But what about the sandwich model, which satisfies 
property A, but not property B? The TBF proof as it 
stands is not immediately valid in these cases. Neverthe- 
less, it turns out that for such models, the TBF proof can 
be made to work by a modification of their structures. 

The basic idea of the modification is as follows. The 
TBF proof of a glassy transition is based on drawing 
structures consisting of sets of overlapping rectangles, 
showing that there is a sufficiently high probability (fi- 
nite for the proof of discontinuity, and approaching 1 for 
the proof of diverging crossover length) that each rect- 
angle has a spanning path in the desired direction, and 
then using property B to conclude that the intersecting 
paths form a frozen cluster. For models that lack prop- 
erty B, we use the same figures as in the TBF proof 
(e.g. Figs 2a and 2b of Ref. [3), but pick the rectangle 
sizes large enough that there is a high probability that 
each rectangle has multiple spanning paths (C(L 1_Z ) for 
a rectangle of O(L)), each occurring in a disjoint parallel 




FIG. 10: The modification in the TBF discontinuity struc- 
ture for models failing the "No Perpendicular Crossing" rule 
(Property B) 

subrectangle. Then in each place where the TBF proof 
assumes a T-junction based on property B, we will have 
a northeast (northwest) path crossing many northwest 
(northeast) paths. The probability that no T-junction 
occurs turns out to decay exponentially with the number 
of northwest (northeast) paths. In figure[lO]we show how 
the discontinuity structure of the TBF proof (from Fig. 
2a of Ref. [14]) is modified by this procedure. 

To implement these ideas, we need to modify Propo- 
sition 5.1 of Ref. [15{ . which says that sufficiently large 
rectangles of size aL x L are very likely at the critical 
point to have chains connecting the sides of length aL: 

Proposition 5.1 (From Ref. [IB] ): 
For L — > oo there exists c > s.t. 

aL (~ northeast occupied cluster connecting the 
sides of length aL) 

< exp(-cL 1 ~ z ) 

Here z is the the dynamical exponent, which has been nu- 
merically found to be approximately 0.63 in two dimen- 
sions [43| . The TBF proof of this Proposition assumes a 
conjecture that there is a well-defined z < 1 (Conjecture 
3.1 of Ref. (HI)- We replace their Proposition 5.1 with 
the following proposition, which instead states that we 
are likely to have 0{L X ~ Z ) connecting chains in disjoint 
parallel subrectangles: 

Proposition 100: 

For L — > oo there exists c > and r > s.t. 

p P £ aL {~^ \rL 1 ~ z \ disjoint northeast occupied clusters, 
occuring in disjoint parallel subrectangles, 
connecting the sides of length aL) 

< exp(— cL 1_z ) 



Proof. We divide the box of size L by aL into aL x ~ z par- 
allel disjoint subrectangles, each of size LxL z . Assuming 
the conjecture of anisotropic scaling in directed percola- 
tion (conjecture 3.1 of Ref. [151]). each subrectangle has 
a probability q > of having a path connecting the two 
sizes of length L z , contained within that subrectangle. 
The expected number of crossings is qaL x ~ z , and for any 
r < qa the probability of having less than rL 1 ~ z crossings 
decays exponentially in L l ~ z . □ 
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The TBF proof of discontinuity shows that for certain 
structures of rectangles, there is a nonzero probability 
that each rectangle has a suitable "event," and that, as- 
suming property B, the existence of each event (a rect- 
angle crossing) results in a stable structure at the critical 
point. Now, with our modified proposition, each "event" 
is the presence of multiple crossings (in disjoint paral- 
lel subrectangles) in each rectangle, rather than single 
crossings. 

Without property B, this does not guarantee a frozen 
cluster. However, we see in appendix [A] that when a 
northeast (northwest) path crosses n northwest (north- 
east) paths, in n disjoint subrectangles, the probability of 
not forming a single T-junction decays exponentially in n. 
This is physically obvious, since for large subrectangles, 
the probability of each T-junction in each subrectangle is 
essentially independent; however, since the probabilities 
are not truly independent, more work is needed to make 
this rigorous. The details of the proof are relegated to ap- 
pendix [AJ More generally, the arguments in appendix |A1 
show that we can treat the probabilities of T-junctions 
in different subrectangles as independent, when establish- 
ing an upper bound on the probability. We will use this 
throughout this section to multiply such probabilities as 
if they were independent. 

Assuming property B, the TBF proof shows that the 
probability of having suitable events in each rectangle 
is nonzero at the critical point. We now need to show 
that, even without property B, this results in a finite 
probability of an appropriate set of T-junctions. In the 
TBF discontinuity structure, shown in Fig. 6 of Ref. , 
the rectangles are labelled 72.^' , i > 1. Each rectangle 
1Z\' 2 is twice as large as the rectangles So the 

"suitable events" of Proposition 100 give 71.' at least 
fc2 s ( 1_z ) crossings parallel to its long direction, where k 
is some positive constant. The arguments in appendix [A] 
show that this results in a T-junction with probability 

(1 — (1 — r) k2 ( 1 ), for some positive r and k. 

Starting at the origin, the probability of forming ap- 
propriate T-junctions off to infinity can then be seen to 
be 



(1 _ (1 _ ri) ^ )2 JJ (1 _ (1 _ r2 )^ )2 ^ (1) 

i=1 

for some positive r%, r2, and k. This product converges 
to a positive number, so the transition is proven to be 
discontinuous (subject to assumption of the well-tested 
conjecture of an anisotropic critical exponent in directed 
percolation) . 

The proof of the diverging crossover length can be 
made to avoid the assumption of property B by a simi- 
lar modification of the TBF structures. Again, we begin 
by repeating the TBF structures, with the set of parallel 
rectangles in figure 2a of Ref. In that picture, if 

every rectangle has a spanning path, and the paths all 
intersect, there will be a spanning frozen cluster. TBF 



F 




E 

FIG. 11: Neighboring relations for the pinwheel model. 



consider the case where each rectangle has sides of order 
the directed percolation parallel correlation length, 
They then show that C3 and C4 can be chosen such that if 
the system size is L < exp[c3^(p) 1_z ], the probabil- 
ity that each rectangle is occupied by a spanning cluster 
approaches 1 as L — > 00, p —> p~ . If property B were to 
hold, this would result in T-junctions that would create 
a frozen structure, and show that the crossover length 
diverges as p — > p~ . We no longer have property B; but 
instead, by replacing Proposition 5.1 with Proposition 
100, we can choose the rectangle sizes such that each 
rectangle is occupied by "many" spanning clusters (with 
"many" defined by Proposition 100). Given this, we can 
start at an arbitrary rectangle, and then work our way 
out, looking for T-junctions to create a spanning frozen 
structure. We will only fail to create a frozen structure 
if at some point we reach intersecting rectangles where 
one spanning path in one rectangle crosses many span- 
ning paths in the other rectangle, but without creating a 
T-junction. By the arguments in appendix A, the proba- 
bility of this occurring decays exponentially in cL 1 ^ z for 
some c > 0. So even with 0(L/£||) 2 intersections, the 
probability that we ever fail in this process goes to as 
L — > 00 and p — > p~ , and we are essentially guaranteed 
a frozen structure. This shows that the crossover length 
diverges as we approach the critical point, with the same 
lower bound that TBF found. 



V. PINWHEEL MODEL AND 8-SPIRAL MODEL 

The models we have discussed so far have two possible 
ways in which a site can be stable, and varying neigh- 
boring relations for the culling relations. However, we 
can also define generalizations in which there are three 
or more ways in which a site can be stabilized. For ex- 
ample, for the Pinwheel model, shown in figure 111! the 
condition for a site to be stable is that it (1) have neigh- 
bors in both the sets A and B, or (2) have neighbors in 
both the sets C and D, or (3) have neighbors in both the 
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FIG. 12: Neighboring relations for the 8-spiral model in which 
there are four processes by which a site can be made stable. 
For a site to be stable it must have neighbors in A and B, or 
C and D, or E and F, or G and H. 



sets E and F. This gives a site three possible directions 
for stabilizing chains. Similarly, figure \TZ\ shows a model 
in which there are four possible directions for a site to be 
stable, such that there is an extra "or" : for the sets G 
and H. We denote this model the 8-spiral model. 

Despite the extra ways in which sites can be made 
stable, the TBF proof of a glassy transition is still valid, 
because Property A holds for both the pinwheel model 
and the 8-spiral model. That is, in both of these models if 
two A-B chains (or two C-D chains, or two E-F chains, or 
two G-H chains) cross, they must have sites in common. 
This turns out to be sufficient to show that there will be 
a void that grows forever in the infinite system limit. 

For the sandwich and spiral models, we needed to show 
that a stable site at the corner of a diamond-shaped void 
has to be part of a directed percolation-like chain (DP- 
like chain) of order the size of the void. For the pinwheel 
model, we consider hexagonal voids, and show that stable 
sites at the corners of these voids must be "associated" 
with a long DP-like chain, where "associated" will be de- 
fined by the construction below. Then, just as for the 
sandwich and spiral models, for p < p c , long DP-like 
chains are exponentially suppressed, giving voids a finite 
probability to grow forever, showing that the infinite sys- 
tem is empty for p < p c . 

Consider the site y at the corner of the hexagonal void 
in figure [T31 Look at all A-B chains coming out of y, 
and pick the lowest possible chain; in other words, look 
for successive A neighbors, and if a site has two possible 
A neighbors, pick the lower one. If that chain reaches 
the dashed line in figure [T3J we have a DP-like chain of 
order the size of the void, and are done. Otherwise, this 
A-B chain must terminate either in a C-D chain, or an 
E-F chain. If it terminates in an E-F chain, that chain 
must terminate in a C-D chain, which must then cross 
the original A-B chain. (Note that the E-F chain cannot 
terminate in an A-B chain, since by the "No Parallel 
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FIG. 13: The top of a hexagonal void in the Pinwheel model, 
and examples of the A-B and C-D chains discussed in the 
text. 



Crossing" rule, the new A-B chain would intersect the 
first A-B chain, and contradict our assumption that we 
chose the lowest A-B chain coming out of the site y.) So if 
the A-B chain coming out of y does not reach the dashed 
line, it must either terminate in a C-D chain or cross a 
C-D chain(the latter case is shown in figure IT5]) . Pick 
the lowest of all the C-D chains that cross or intersect 
our A-B chain. This C-D chain must reach the dashed 
line, using the same logic as before (if it terminated in an 
A-B chain, that would intersect the first A-B chain, and 
contradict the assumption that we chose the lowest A-B 
chain coming out of y; while if it terminated in an E-F 
chain, that E-F chain would have to turn into either an A- 
B or C-D chain before reaching the void, again resulting 
in a contradiction.) Since the A-B and C-D chains that 
we have constructed cross, and together include both y 
and the dashed line, at least one of the chains must be of 
order the size of the void. A similar construction can also 
be used for the 8-spiral model using a diamond-shaped 
void. 

Having established where the critical point is, we fol- 
low TBF, and as before consider an infinite sequence of 
two types of rectangular regions, IZ1.2 (see Fig. ITU)) , one 
of which contains AB paths, and the other of which con- 
tains EF paths. (Any pair of types of paths — AB/CD, 
AB/EF, or CD/EF — are permissible, as are any of the 
6 possible pairs for the 8-spiral model.) Again, the se- 
quence is constructed such that the AB paths and the 
EF paths are mutually intersecting with a frozen back- 
bone that contains the origin, and the TBF proof simply 
carries through with the additional modification we have 
introduced for models that do not obey Property B. All 
in all, since having more neighboring relations gives more 
ways for an occupied site to be stable (without depressing 
the DP critical point), the TBF constructions of discon- 
tinuous percolation structures at the critical point sim- 
ply carry through. Just as for the sandwich model, the 
TBF proof needs to be modified to deal with the lack of 
Property B. One needs at least two intersecting rectangu- 
lar regions in which the probability for two "transverse" 
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blocking directions each undergoing a directed percola- 
tion transition independently is nonzero. 



VI. MODELS IN HIGHER DIMENSIONS 

Consideration of the "No Parallel Crossing" rule shows 
that for higher-dimensional generalizations of the knights 
model, the TBF proof cannot be generalized in a straight- 
forward manner to show provably glassy transitions. We 
will show that the "No Parallel Crossing" rule never holds, 
so the critical point is always depressed below that of di- 
rected percolation. 

To be specific, suppose that in three dimensions we 
have disjoint neighboring sets A, B, C, and D, general- 
izing the northeast, southwest, northwest, and southeast 
sets of the knights model (see figure [T]) . Each set should 
consist of three linearly independent vectors, and the sets 
A and B should be opposite of each other, as should the 
sets C and D. Then, just as in the knights model, the 
culling rule is that a stable site should have occupied 
neighbors in both the sets A and B, or in both the sets 
C and D. Then, just as for the sandwich and spiral 
models, there are two directed percolation processes by 
which a site can be made stable (A — B chains and C — D 
chains), and one might think that for an appropriate set 
of neighboring relations the TBF proof could be used to 
show a glassy transition. 

However, it turns out that because these models never 
satisfy three-dimensional generalizations of the "No Par- 
allel Crossing" rule, the critical point is always depressed 
below that of three-dimensional directed percolation, and 
the TBF proof cannot be directly generalized for these 
models. 

Property A says that two chains running in similar di- 
rections cannot cross without having sites in common. 
For certain two-dimensional models, such as the sand- 
wich and spiral models, these property can be required 
by the topology and neighboring relations. However, in 
three dimensions, the topology always makes it easy for 
two directed chains to miss each other, and so no three- 
dimensional generalization of property A can be satisfied, 
regardless of the neighboring relations. Furthermore, if 
the two chains miss each other, then the buttressing of 
each type of chain does not occur and the resulting tran- 
sition may be continuous. 

This rough argument can be formalized by show- 
ing that three-dimensional generalizations of the knights 
model always have substructures such as the one shown 
in figure[5J That is, it is always possible to find substruc- 
tures that have no long A — B chains connecting their 
ends, but which can join up two A — B chains and sta- 
bilize their ends. Then, by the arguments in section [TTTT 
we should expect these to depress the critical point below 
that of three-dimensional directed percolation. 

Specializing to three dimensions for convenience, let 
A = {0*1,0*2,03} consist of three linearly independent 
3- vectors, and C = {ci, 02,03} consist of three lin- 



early independent 3- vectors with A / C. Also, define 
B = {— a*i, —02, —0*3} and D = {— ci, — c*2, — c 3 }. For an 
occupied site r to be stable, it must have either (1) occu- 
pied neighbors from both r + A and r + B or (2) occupied 
neighbors from both r + C and r + D. The first condi- 
tion we denote the A-B condition, the second, the C-D 
condition. 

If we only enforced the A-B condition, then we would 
just have three dimensional directed percolation (modulo 
finite clusters). However, for the model defined above, 
there is an extra way to be stable, resulting in p c < p® p ' . 
Consider a finite structure with the following properties: 
(1) all occupied sites are stable under the culling rules 
except occupied sites ?i and r/, (2) r% has an occupied 
neighbor in r + B and ff has an occupied neighbor in 
r + A, and (3) there is no AB path connecting r*; and f/, 
but the structure is stable because there exists a path 
where at least one occupied site is stable under the C-D 
condition. We relegate to appendix B the proof of the 
existence of such a finite structure. 

With the finite structures defined above, some occu- 
pied sites that were unstable under the A-B condition 
now become stable. Slightly below p® p , the system is 
about to percolate using the A-B condition alone. The 
substructures act as extra, local bonds, joining up long 
A-B paths and pushing the system above the critical 
point, just as in the knight model. Again, arguments 
similar to these have been made rigorous by Aizenmann 
and Lebowitz in the case of undirected percolation, and 
perhaps can be made rigorous in the case of directed (ori- 
ented) percolation. 



VII. DISCUSSION 

The discovery of a two-dimensional percolation tran- 
sition, where the sudden emergence of a discontinuous 
backbone coincides with a crossover length diverging 
faster than a power law, is recent, and of great interest 
for glassy systems, jamming systems, and phase transi- 
tions in general. Unusual transitions have been found 
previously in mean field systems of a slightly different 
nature, but not in finite dimensions. And while the finite- 
dimensional transition is discontinuous, it is not driven 
by nucleation, as with ordinary discontinuous transitions, 
but instead by a scaffolding of many tenuous directed per- 
colation paths occurring simultaneously to form a bulky 
structure. We have shown that Property A is required 
for the proof that p c = p,P p , but Property B is not. All 
that is needed to prove that the transition is discontinu- 
ous (once p c = p^ p is established) is a finite probability 
for two transverse percolating structures to intersect, to 
prevent each other from being culled. Therefore, one can 
construct other models, such as the pinwhcel and 8-spiral 
models, that exhibit a similar transition in two dimen- 
sions. The phenomena is not as specific as might seem 
at first glance. However, such a buttressing mechanism 
in dimensions higher than two is more difficult because 
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it is more difficult for percolating paths to intersect and 
form a buttressing, bulky structure. 

Models like the knights model, where p c is most likely 
less than p® p , provide physicists, mathematicians, and 
computer scientists with a motivation to study new mod- 
els of correlated percolation — models that are not isomor- 
phic to directed percolation, but quite possibly in the 
same universality class. Once this avenue is pursued fur- 
ther, one can then easily extend the class of models for 
which a finite-temperature transition can be rigorously 
shown. To begin, it would be interesting to consider a 
directed percolation model in two dimensions where the 
number of nearest neighbors is greater than two. For ex- 
ample, if the number of nearest neighbors was increased 
to four, would the percolation transition still be in the 
same universality class as directed percolation? If so, as is 
presumably the case, then one could construct a jamming 
percolation model with sets larger than two sites. These 
jamming percolation models would then be isomorphic to 
the next-neighbor directed percolation models and then 
one could use results from directed percolation to prove 
a percolation transition. 

There exists another class of correlated percolation 
models called force-balance percolation. The first model 
in this class was defined in Ref. [37[ . Other force-balance 
percolation models are currently being constructed and 
studied (38|. The force-balance percolation models dif- 
fer from the jamming percolation models in that (1) the 
sets, such as A, B, etc., are overlapping and (2) the 
"or" between pairs of sets is changed to "and". Given 
these differences, the methods of proof used here cannot 
be easily applied. Numerical results indicate that the 
transition is discontinuous with a nontrivial "correlation 
length" exponent, indicating that the transition may not 
the garden-variety discontinuous transition. The force- 
balance models are perhaps less artificial in that they 
mimic force-balance by requiring that an occupied site 
(i.e. a particle) have occupied neighbors to its left and 
right, as well as its top and bottom, in order to be stable. 
However, little has been rigorously proven about them. 
To make progress along these lines would be useful. 

What about the lack of finite stable clusters in models 
of jamming percolation? Recent numerical work on an- 
other correlated percolation model, fc-core percolation, 
with k = 4 on the four-dimensional hypercubic lattice, 
appears to exhibit an ordinary, discontinuous percola- 
tion transition driven by nucleation [44| . Finite clusters 
exist in this model, unlike in the jamming percolation 
models. Therefore, in the jamming percolation models, 
there can be no surface tension between the percolating 
and nonpercolating phases, which is typical of an ordi- 
nary discontinuous transition. If finite clusters are al- 
lowed in a correlated percolation model, one might guess 
that the unusual nature of the transition would be de- 
stroyed. However, finite clusters, other than individually 
floating particles, do not appear in the jamming transi- 
tion of granular particles. Otherwise, the packing would 
not be static. So it is unclear whether the existence of 



finite clusters pertains to the jamming of granular parti- 
cles. This is also the case for the glass transition. 

If more of an analogy between jamming and models of 
jamming percolation is to be made in finite dimensions 
(setting aside the matter of the critical dimension of jam- 
ming), a model where the fraction of sites participating 
in the infinite cluster increases smaller than linearly just 
above the transition must be found. Furthermore, the 
existence of a jamming percolation model with a univer- 
sal jump in the number of occupied sites at the transition 
that "naturally" emerges as opposed to being externally 
imposed [451 ] . is yet another necessary quest if a jamming 
percolation model of jamming is to be found. 

Finally, models of correlated percolation, such as the 
sandwich model, tell us that there do indeed exist kincti- 
cally constrained models of glassy dynamics that exhibit 
unusual phase transitions in finite dimensions. Therefore, 
this avenue of exploration for understanding possible 
finite-temperature glass transitions in finite-dimensions 
remains open. Since our work helps to clarify which 
jamming/correlated percolation models can be rigorously 
shown to have an unusual finite-temperature glass tran- 
sition with a particular set of properties, other models 
exhibiting possibly other unusual behaviours can hope- 
fully be more easily developed in the near future. 



APPENDIX A: PROOF THAT MULTIPLE 
CROSSINGS ARE EXPONENTIALLY UNLIKELY 
TO AVOID CREATING T-JUNCTIONS 

In this appendix we justify the claim made in sec- 
tion IIVI that given n crossings, the probability that no 
crossing results in a T-junction decays at least expo- 
nentially in n. This would be immediately true if each 
crossing resulted in an independent probability of a T- 
junction. So what we show is that these crossings, by 
occurring in disjoint subrectangles, can be effectively 
treated as independent (in establishing an upper bound 
on the probability). 

The relevant picture is shown for n = 4 in Fig. 1141 
There are n disjoint rectangles, labeled by i, 1 < i < 
n, each of which has at least one northwest spanning 
path. We call this event H, and conditionalize upon the 
occurrence of H . Each northwest path must cross the 
northeast spanning path, but without Property B, these 
crossings do not necessarily result in T-junctions, where 
by a T-junction we mean specifically a site in common 
between the paths that stabilizes the northeast part of 
the northeast path. Let Gj be the event that at least one 
crossing in rectangle i forms a T-junction. 

Now for a configuration with ->Gj, look at the sites 
in the vicinity of the crossing (if there is more than one 
crossing, we choose one by an arbitrary ordering of possi- 
ble crossing locations). Restricting ourselves first to the 
sandwich model, Fig. 1151 shows the only way that Gj can 
fail to happen. The site labeled by a star must be vacant, 
and if that site is made occupied, the new configuration 
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FIG. 14: Multiple crossings, resulting in at least one In- 
junction with high probability. Crossings are drawn with 
open circles to emphasize that these crossings may or may 
not result in sites common to the crossings paths. 



/i'(-.Gi | H) < c^(G l | H) 



(A2) 



for some < c < 1. 

Since the configuration changes only take place within 
a rectangle i, they do not affect whether or not we have 
Gj for j ^ i, and we can write the above inequality for 
a state where we specify whether these other Gj occur. 
For example, for i — 3 we might write 



m"(Gi n -.g 2 n -.g 3 n G 4 • ■ ■ n G„ | H) < 

c^ p (Gi n -G 2 n G 3 n G 4 • • • n G n \H) , (A3) 

with the same c as above. This can be intuitively thought 
of as treating the different probabilities of forming T- 
junctions as independent. Repeatedly using equation IA3I 
in different subrectangles, we find that 
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FIG. 15: The only configuration in the sandwich model where 
northeast and northwest chains cross, but the northeast part 
of the northeast chain is not stabilized by a T-junction. If 
the site marked with a star was occupied, this would create a 
T-junction. 



is in Gi. So a local change in the vicinity of the cross- 
ing can always create a T-junction. This local change 
induces a mapping /, from the set of states with ^Gj to 
a subset of the set of states with G;. The mapping is one- 
to-one onto this subset, and for any state S with ->Gj, the 
probability of the configuration fi(S) is p/(l — p) times 
the probability of the configuration S, where p is the site 
occupation probability. Thus 



/AG I H) 



(Al) 



More generally, if we want to consider other variations 
of the knights model that lack Property B, we need only 
that for any crossing without a T-junction, some local 
configuration of changes in a bounded region around the 
crossing can create a T-junction. The induced mapping 
can be many-to-one, so long as the "many" is bounded 
(which follows automatically from the restriction that the 
configuration changes occur in a bounded region around 
the crossing). This will more generally give 



"(.Or G,|ff ) £ ar: 



(A4) 



We are thus exponentially unlikely to have no T- 
j unctions. 



APPENDIX B: EXISTENCE OF FINITE 
STRUCTURES IN THREE-DIMENSIONAL 
MODELS 

In this appendix we prove the existence of the finite 
structures discussed in the section on jamming percola- 
tion in three-dimensions. Before giving the formal proof, 
we sketch the qualitative idea behind the construction of 
these structures. In the structure for the knights model, 
shown in figureO there are two parallelograms consisting 
of A — B (northeast-southwest) chains. The two parallel- 
ogragrams are parallel to each other, and are connected 
and stabilized by G — D (northwest-southeast) chains. 
In two dimensions, such a figure can only be constructed 
by having some chains cross each other, and this results 
in a long A — B chain connecting the two ends, unless 
the model violates Property A. However, in three dimen- 
sions, the C — D chains can always be run in a direction 
independent of the plane of the A — B parallelograms, so 
such a substructure can always be formed, regardless of 
the neighboring relations. 

We now formalize this argument. Recalling that a±, a 2 , 
S3, and ci are four three-dimensional vectors, with the 
first three being linearly independent and c\ not equal to 
any of the first three, there must exist mi, m 2 , 7713 , and 
n such that 
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FIG. 16: The structure in the case where ci = 2i=i m *<^ 
with all mi > 0. 



Since the Si are linearly independent of each other, n =/= 0. 
Since c\ A, at least two of the m, must be nonzero. 

If all of the rrii are positive, then it is easy to make 
the desired structure. Let F = Xh=i m i®i — n ^i- Then 
make the structure in figure [TBI where each of the vectors 
F represents an A-B chain of length X)i=i m ii ano - nc i 
represents a C-D chain of length n. This structure has 
the desired properties. If all of the rrii are negative, we 
simply replace C with D and use the same argument. 

We are now left with the most complicated case, where 
some rrii are positive, and some rrii are negative. In this 
case, we redefine the rrii, and rewrite equation [BTJ as 

3 3 

^ m t a t = ™i<ii + "Ci ! ( B2 ) 

i=l i=l 

where n, all rrii, and all rhi, are positive, and for any i, 
either rrii or m, is zero. Then define F = X)i=i an( i 
G = X)i=i WH^i" F and G are both nonzero and linearly 
independent. 

We can now make the structure shown in figure [T7] 
The vectors 2F and 2G represent A-B chains of length 
2^i=i m i an d 22 j=1 mj. In these chains every site is 
stable by the A-B condition except for fi and r\. The 
vectors nc\ are chains of length n in which every site is 
stable by the C-D condition. 

Next, since {ai, 02, 03} form a complete basis for three- 
dimensional space, and {F, G} are only two vectors, there 
exists a vector H = ^ =1 p»Oi, with all pi > 0, such 
that H is linearly independent of {F, G}. (Note: this 
is where the three-dimensional case differs from the two- 
dimensional case. For models such as the spiral and sand- 
wich models, the two vectors F and G already span the 
space.) We can then make a second copy of figure [TTJ 
displaced from the original by H , as shown in figure [T51 

We can now check that in figure [TS] all sites except 
the start and end sites, fi and ry are stable under the 
culling rules, and that fi and ff have neighbors from B 
and A, respectively. It remains to check that there is no 
A-B chain connecting f to ff in this structure. There 
is no obvious such chain, but depending on the vectors 




FIG. 17: A substructure 




FIG. 18: The full structure with the desired properties 



Si and ci, it is possible that there are some sites in the 
F, G, and H chains by chance are separated by a vector 
Si, inadvertently forming an A-B chain between r*, and 
ff. However, if this is the case, we can create new larger 
structure, simply by multiplying n, all rrii, all f^i, and 
all f>i by the same multiplicative constant. The structure 
thus grows larger, while the vectors stay the same, 
so for a sufficiently large multiplicative constant, it is 
impossible for the different chains in the structure to be 
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adjacent by Si connections. We thus have a structure with the desired properties. 
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